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Abstract 

I study the critical behavior of a two-dimensional dimer-trimer lattice model, 
introduced by Kohler and ben-Avraham [J. Phys. A 24, L621 (1991)], for het- 
erogeneous catalysis of the reaction + \Bz —* AB. The model possesses 
infinitely many absorbing states in which the lattice is saturated by adsorbed 
particles and reactions cease because only isolated vacancies are left. Results 
for various critical exponents show that the model exhibits the same critical 
behavior as directed percolation, contrary to earlier findings by Kohler and ben- 
Avraham. Together with several other studies, reviewed briefly in this article, 
this confirms that directed percolation is the generic universality class for models 
with infinitely many absorbing states. 



PACS Numbers: 05.70.Ln, 05.50. +q, 64.90. +b 



1 Introduction 



The critical behavior of nonequilibrium lattice models with absorbing states has at- 
tracted a great deal of interest in recent years. It has been shown that most models 
exhibiting a continuous phase transition to a unique absorbing state belong to the same 
universality class. The best known examples are probably directed percolation (DP) 
U, Reggeon field theory 0, [|, the contact process ||, |J, and Schlogl's first and second 
models |§— 1§. Extensive studies of these and many other models [f|-[|16| have revealed 



that the critical exponents are unaffected by a wide range of changes in the evolution 
rules. This provides firm support for the DP conjecture, first stated by Janssen |7| 
and Grassberger || , that directed pecolation is the generic critical behavior of models 
exhibiting a continuous transition into a unique absorbing state. 

The study of models with more than one absorbing state is in an early stage, and there 
is still some controversy regarding the possible universality classes for such models. 
Models with infinitely many absorbing states arise naturally in the study of reactions 
catalysed by a surface as soon as the absorption mechanism for the various species 



requires more than one vacant site. Some studies of such models |17], [18], [T{J in two 
dimensions indicated that they did not belong to the universality class of directed per- 
colation. However, the exponent estimates varied significantly from model to model 
and could therefore not be consistent with just a single new universality class. While 
the exponent estimates quoted in some studies clearly ruled out DP critical behavior, 
those of other studies differed only marginally from the DP values. On the other hand, 
studies of several one-dimensional models |2D], [H], clearly placed these in the DP 
universality class, at least as far as the static critical behavior is concerned. In all of 
these models the number of absorbing configurations grows exponentially with system 
size, though the absorbing configurations are characterized by the vanishing of a unique 



quantity, e.g., the number of particle pairs |2(J or in other cases [[17], [L8| the number 
of nearest neighbor vacancy pairs. Although there still is some evidence suggesting 
that models with infinitely many absorbing state may exhibit non-DP critical behav- 
ior, generally it seems that such models belong to the DP universality class. In this 
article I present results from a study of the dimer-trimer model, introduced by Kohler 



and ben-Avraham [T7[, which shows that this model, despite earlier evidence to the 



contrary, also belongs to the DP universality class. This result lends further support to 
the extensions of the DP conjecture to models with multiple components [23] and/or 



infinitely many absorbing states |20], |21| , at least in cases where the absorbing states 



are characterized by the vanishing of a unique quantity. 

The remainder of this article is organised as follows. In Section I will briefly present 
the various models and review the results of previous studies. I will present the results 
from my study of the dimer-trimer in Section |3] and summarise and discuss the evidence 
regarding universality in Section f|. 
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2 Models with infinitely many absorbing states. 



The dimer-trimer lattice model for heterogeneous catalysis for the reaction -^K^ + 
IB3 — > AB was introduced by Kohler and ben-Avraham |P7|| . Adsorption of dimers 
(trimers) is attempted with probability p (1 — p) and succeeds if the molecule hits a 
pair (triplet) of nearest neighbor empty sites. Immediately upon adsorption the dimer 
(trimer) dissociates and each site of the pair (triplet) becomes occupied by one A (B) 
particle. If A and B particles happen to become nearest neighbors they react and 
the AB product desorbs at once, leaving behind two empty sites. Configurations with 
only isolated empty sites and adsorbed particles are absorbing; the number of such 



configurations grows exponentially with lattice size. Simulations [|I7| revealed that 
when p < pi the system enters a trimer-saturated state with only adsorbed B particles 
and isolated empty sites. Likewise a dimer-saturated state is reached for p > p 2 , and 
only for intermediate values does the system possess an active steady state in which the 
production of AB goes on indefinitely. The transition at p\ = 0.3403(2) is continuous 
whereas the transition at p 2 = 0.4610(8) is discontinuous. Critical exponents (3a and 
(3b describing, respectively, the behavior of the density of A and B particles near p 1 may 
be defined as: pa oc \p — p\\^ A and p^ 1 — Pb oc \p — Pil 13 * 3 , where p^ is the saturation 
concentration of B at p\. Steady state computer simulations |DJ yielded (3a = 0.80(6) 
and (3b = 0.63(5). For directed percolation in (2+l)-dimensions (3 = 0.592(10) (this 
estimate is obtained using the scaling relation (3 = 5^ || with 5 = 0.460(6) and 
v\\ = 1.286(5) [p5|0. While the estimate for (3a is well above the DP value the estimate 
for (3b is marginally consistent with it. The dimer-trimer model was also studied using 
time-dependent simulations. The general idea of time- dependent simulations is to 
start from a configuration which is very close to the absorbing state, and then follow 
the "average" time evolution of this configuration by simulating a large ensemble of 
independent realisations 0. Out of the many possibile (near-absorbing) initial states, 
Kohler and ben-Avraham used absorbing states generated by the critical system, and 
then placed a triplet of vacancies at the center (see the next section for more details). 
As usual in this type of simulation they measured the survival probability P(t), the 
average number of empty sites n(t), and the average mean square distance of spreading 
R 2 (t) from the origin. These quantities are expected to exhibit power-law behavior 
at p\ characterized by exponents —5, i] and z, respectively. The estimates for the 
dimer-trimer model fll7 |, 5 = 0.40(1), rj = 0.28(1), and z = 1.19(1), clearly differ from 



the DP values 5 = 0.460(6), r] = 0.214(8) and z = 1.134(4). These results led 
to the conclusion that the dimer-trimer model belongs to a new universality class. As 
already mentioned the results of my study of this model (reported in Section |3p do 
not support this conclusion but shows that the dimer-trimer model belongs to the DP 
universality class. 



The dimer-dimer (DD) model [|T8j is based on the oxidation of hydrogen on a metal 
surface. 2 adsorption is attempted with probability p, and H 2 adsorption with proba- 
bility 1— p. Both O2 and H 2 require a nearest neighbor pair of vacancies; both dissociate 
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upon adsorption. Whenever H and are nearest neighbors they react to form OH, 
which resides at a single site. Similarly, OH reacts with neighboring H, forming H 2 
which desorbs immediately. In addition adsorbed H atoms are allowed to diffuse on the 
surface. The DD model has been studied in several versions including further processes 
such as, recombination and desorption of H 2 , desorption of OH, and reactions between 
neighboring OH molecules (leading to the formation and desorption of H 2 leaving 
behind one adsorbed O atom). Depending on which processes are included the DD 
model exhibits one or two continuous phase transitions. For p < pi, the steady state 
is absorbing and comprised of a mixture of O, OH and isolated vacancies. In certain 
versions of the model the lattice becomes saturated with H for p > p 2 - a unique 
absorbing state. For p x < p < p 2 there is an active steady state with ongoing produc- 
tion of H 2 0. Monte Carlo simulations |1| generally yielded estimates for flx — 1/2 
with an uncertainty of approximately 5-10%, where flx describe the power-law decay 
of various densities at p\. However, in certain cases it was found |18j that (3x — 2/3 or 
1. Most of the results are not consistent with directed percolation and could indicate 
that the DD model belongs to a new universality class. However, one would expect the 
exponents, corresponding to different densities, to be equal. The difference in some of 
the reported values thus suggests substantial statistical uncertainties. 



Jensen and Dickman studied two simpler models, the pair contact process (PCP) and 
the dimer reaction (DR) model, presenting infinitely many absorbing states in the 
simpler context of single-component, one-dimensional models PU| , |2T[ . In the PCP, 
nearest neighbor pairs of particles annihilate mutually (with probability p) or else 
(with probability 1 — p) create a new particle at an empty nearest neighbor. Any 
configuration without pairs is absorbing; there are evidently many such states (> 2 N / 2 
on a 1-d lattice of N sites). Steady state simulations, including a thorough finite- 
size scaling analysis, yielded exponent estimates consistent with DP critical behavior. 
In the one-dimensional DR, particles are not allowed to occupy neighboring sites. If 
sites i, i — 1 and i + 1 are vacant adsorption may take place at site i. Suppose a 
particle has just arrived at site i. If sites i — 3, i — 2, i + 2, and i + 3 are all vacant, 
the particle remains. If any of the four sites are occupied, the new particle reacts 
with one other particle with probability 1 — p, and remains with probability p. The 
second neighbors have priority in the reaction. Any configuration without triplets of 
empty sites is absorbing. Again steady state simulations (including finite-size scaling) 
revealed critical behavior in the DP universality class. 

Another single-component model which exhibits a continuous phase transition from an 
active steady state to one of a multitude of absorbing configurations is the threshold 
transfer process (TTP) |[22|| . In the TTP sites may be vacant, or singly or doubly 
occupied, which can be described by an occupation variable <7j = 0, 1 or 2. If o"j = 
0,then o~i — > 1 with probability p\ if Oi = 1, then o~i —>■ with probability 1 — p. In 
the absence of doubly-occupied sites, the dynamics trivially leads to a steady state 
in which a fraction p of the sites have o~i — 1. When o~i = 2 particles may move to 
neighboring sites. If <Tj_i < 2, one particle moves to that site; likewise, a particle 
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moves from i to i + 1 if o~j + i < 2. cr, is diminished accordingly in this deterministic, 
particle-conserving transfer. Configurations devoid of doubly-occupied sites form an 
absorbing subspace with trivial dynamics and can be avoided only if p is sufficiently 
large. Simulations |22[ revealed that the steady state critical behavior placed the TTP 



in the DP universality class. 

Recently, Yaldram et al, studied the critical behavior of a model for the CO+NO— > 
C0 2 +|N 2 catalytic surface reaction. With probability p a CO molecule is adsorbed on 
an empty site and with probability 1 — p NO adsorption is attempted. NO dissociates 
upon adsorption and therefore requires a nearest neighbor pair of empty sites. After 
each adsorption the nearest neighbors are checked (in random order) and CO+O reacts 
to form CO2 which leaves the surface at once, likewise N+N forms N 2 which desorbs 



immediately. Computer simulations by Yaldram et al. |19fl showed that when p < pi 
the system enters an absorbing state in which the lattice is covered by a mixture of 
O and N. Again, the number of absorbing configurations grows exponentially with 
system size. At p\ the model exhibits a continuous phase transition into an active 
state in which the catalytic process can prodeed indefinetely. Finally when p exceeds a 
second critical value P2 the model exhibits a discontinuous phase transition into a CO 
and N covered state. Near the critical point p\ one expects the concentrations px of 
various lattice sites X (X = O, N, CO, or an empty site) to follow simple power laws, 
px — p s x {p~Pi) l3x i where p s x is the saturation concentration. Yaldram et al. [|K 



found that j3x = 0.20 — 0.22, which is much smaller than the DP value. This could 
indicate that the CO-NO model belongs to a new universality class. However a more 
thorough study [pj| yielded exponent estimates consistent with DP critical behavior. 



3 Results for the dimer-trimer model. 



In this section I report the results of extensive time-dependent and steady state sim- 
ulations of the dimer-trimer model including a thorough finite-size scaling analysis. 
The model was studied on a triangular lattice with each site either empty or occupied 
by a single A or B particle. In the actual simulations dimer (trimer) adsorption is 
attempted with probability p {1—p). Close to the critical value p\ the lattice soon be- 
comes covered with B-particles. Since all processes depend on the presence of nearest 
neighbor pairs or triplets of empty sites an efficient algorithm uses a list of non-isolated 
empty sites or 'active' sites. By chosing the first site (randomly) from the list of active 
sites the time spent on failed adsorption attempts is greatly reduced. Dimer adsorp- 
tion proceeds via chosing (at random) a nearest neighbor of the active site and placing 
A-particles on these sites if the second site is empty. For trimer adsorption two nearest 
neighbors of the actice site are chosen such that the sites constitute the vertices of an 
equilateral triangle and adsorption takes place if all sites are empty. After each at- 
tempted adsorption the time variable is incremented by 1/N a , where N a is the number 
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of active sites prior to the attempt. Each time step thus equals (on the average) one 
attempted update per lattice site. After each successful adsorption the neighbors of 
newly adsorbed particles are checked to see if A-B pairs were formed and such pairs 
are removed. If an A or B particle is part of more than one pair a random choice is 
made between the different A-B pairs. 



3.1 Time-dependent behavior. 

Time-dependent simulations || is generally the most efficient way of obtaining a pre- 
cise estimate for the location of the critical point in models with absorbing states. The 
general idea is to study the average evolution of a system which initially is very close 
to an absorbing state. For models with a unique absorbing state this is a very simple 
procedure yielding accurate estimates for both the critical point and various critical 
exponents describing the critical power-law behavior of quantities such as the survival 
probability or the average number of active sites. For models with multiple absorb- 
ing state the situation is more intricate. A recent study |21| showed that the critical 



exponents depend upon the choice of initial configuration. However, two important 
facts emerged from this study, first of all the value of the critical point was always 
predicted correctly, and secondly by using an initial configuration generated by run- 
ning the system at the critical point starting from an empty lattice the predictions for 
the dynamical critical exponents coincide with those expected from the static critical 
behavior. A recent more thorough study by Mendes et al. [^2| confirmed this picture 
and led to a generalized scaling ansatz for models with multiple absorbing states. In 
this study I generate the initial configuration by simulating the dimer-trimer model 
on a 128x128 lattice (with periodic boundary conditions) at the value of p under 
investigation until it enters an absorbing state. An off-set (x, y) is then chosen ran- 
domly on this lattice. Hereafter the configuration is mapped cyclically onto a larger 
(512x512) lattice such that the state of site on the large lattice is the same as 
that of site (i + x mod 128, j + y mod 128) on the small lattice. Hereafter a triplet 
of nearest neighbor empty sites is placed at the origin. The size of the large lattice 
ensures that the cluster of empty sites grown from the seed at the origin never reaches 
the boundaries of the lattice. We thus start in a configuration close to an absorbing 
state (just three sites are open) and it should be close to a typical absorbing state 
of the infinite system. For each value of p I simulated 50 independent configurations 
and for each such configuration I simulated 5000 independent samples for a total of 
250,000 samples. Each run had a maximal duration of 2000 time steps, but most sam- 
ples entered an absorbing state before this limit was reached. I measured the survival 
probability P(t) (the probability that the system has not entered an absorbing state at 
time t), the average number of active sites n(t), and the average mean square distance 
of spreading R 2 (t) from the center of the lattice. Notice that n(t) is averaged over all 
runs whereas R 2 {t) is averaged only over the surviving runs. In accordance with the 
scaling ansatz for models with a unique absorbing state || |25|] one expects that these 
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quantities have the following scaling form, 



P(t) oc r*$(At 1 /"n) j (l) 
n(t) oc ^*(At 1/l/ H), (2) 
i? 2 (t) oc t 2 6(At 1/ly ii), (3) 



where A = \p — p±\ is the distance from the critical point, and v\\ is the correlation 
length exponent in the time direction. If $, and are non-singular at the origin 
then asymptotically (t — > oo) Pit), n(t), and R 2 (t) behave as power-laws at pi with 
critical exponents —8, rj, and z, respectively. Generally there are corrections to the 



pure power law behavior so that, e.g., P(t) is more accurately given by |25| 



P(t) oc r 5 (l + ar l + br s ' + ■■■) (4) 

and similarly for fi(t) and R 2 (t). More precise estimates for the critical exponents can 
be obtained if one looks at local slopes 



log(m) 

and similarly for rj(t) and z(t). In a plot of the local slopes vs 1/t the critical exponents 
are given by the intercept of the curve for p± with the y-axis. The off-critical curves 
often have very notable curvature, i.e., one will see the curves for p < p\ veering 
downward while the curves for p > p\ veer upward. This enables one to obtain accurate 
estimates for p\ and the critical exponents. In Fig. 1 I have plotted the local slopes 
for various values of p. From the plot of r](t) it is clear that the two lower curves, 
corresponding to p = 0.3418, and 0.3420, veer downward showing that p\ > 0.3420. 
Likewise the upper curve, p = 0.3424, has a pronounced upward curvature. I therefore 
conclude that p\ = 0.3422(2). This estimate differs quite a bit from that of Kohler and 
ben-Avraham (pi = 0.3403(3)), which is probably due to slightly different algorithms. 
From the intercept of the critical curves with the y-axis I estimate 5 = 0.46(1), rj = 
0.225(5) and z = 1.13(1). These values agree very well with those obtained from 



computer simulations of directed percolation in (2+l)-dimensions [25|, 5 = 0.460(6), 
r] = 0.214(8) and z = 1.134(4). 

From these results it seems reasonable to conclude that the dimer-trimer model belongs 
to the DP universality class. However, due to the somewhat arbitrary choice of the 
initial configuration employed in the time-dependent simulations it would be nice to 
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validate this conclusion through other means. To this end I have also performed 
extensive steady state simulations using a finite-size scaling analysis. 



3.2 Finite-size scaling behavior. 

Finite-size scaling, though originally developed for equilibrium systems, is also applica- 
ble to nonequilibrium second-order phase transitions as demonstrated by Aukrust et al. 
TI|| . As in equilibrium second-order phase transitions one assumes that the (infinite- 
size) system features a length scale which diverges at criticality as, £(p) oc A~ u± , 
where v± is the correlation length exponent in the space direction. The basic finite- 
size scaling ansatz is that the various quantities depend on system-size only through 
the scaled length L/£, or equivalently through the variable AL 1 ^ U± , where L is the 
linear extension of the system. Near the critical point p\ one would expect the steady 
state concentrations px of various lattice sites X, X = A, B, active sites or empty 
sites (note that empty sites include only the isolated vacancies) , to follow simple power 
laws, 



\pf - p x \ oc (p - Pl f\ (6) 

where p s x is the saturation concentration. Note that the saturation concentration for 
active sites and A is zero, whereas it is non-zero for B and empty sites. Thus we 
assume that the density of various sites depends on system size and distance from the 
critical point as: 



\pf(p,L) -p x (p,L)\ oc L-^^iAL 1 ^), (7) 
such that at the critical point pi 

\p s /(px,L)-p x (p 1 ,L)\^L-^\ (8) 

Px, and other quantities, are averaged over the surviving samples only. Fig. 2 shows 
a plot of the average concentration of sites log 2 \p x t {pi,L) — px(pi,L)\ as a function 
of log 2 L at the critical point, p\ = 0.3422. All simulations were performed on lattices 
of size L x L using periodic boundary conditions. The maximal number of timesteps 
in each trial, £m, and number independent samples, N s , varied from tu = 150, N s = 
50, 000 for L = 8 to t M = 75, 000, N s = 500 for L = 256. The slope of the line drawn 
in the figure is f3/v± = 0.81, which comes from the DP estimate (3/v±_ = 0.81(2), using 



the earlier cited estimate for (3 and u± = 0.729(8) [25||. The data falls very nicely on 
the line drawn using the DP estimate thus confirming that the model belongs to the 
DP universality class. 
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Near the critical point the order parameter fluctuations grow like a power law, x = 
L d ((p 2 ) — (p) 2 ) oc A 7 , from which we expect the following finite-size scaling form, 



X (p,L)^L^g(AL 1 ^), (9) 

such that at pi 

X (pi,L)ocL^. (10) 

Fig. 3 shows a plot of log 2 [x(pi, L)] vs log 2 L, as obtained from the fluctuations of the 
number of active sites. The slope of the straight line is 0.39 as obtained from the DP 
value 7/Vj_ = 0.39(2), where I used that 7 = 7 DP -z/|| = 0.285(11) with 7 DP = 1.571(6) 



25]. The excellent agreement between the data and the DP-expectation confirms the 



DP critical behavior of this model. 

One expects a characteristic time for the system, say the relaxation time, to scale like 

r(p, L) oc L-V^TCAL 1 /^), (11) 

such that at pi 



T 



(p^ocL-V-x. ( i2) 



In Fig. 4 I have plotted log 2 [r/ l (pi, L)], where is the time it takes for half the samples 
to enter an absorbing state, as a function of log 2 L. The slope of the line drawn in 
the figure is v\\/v±_ = 1.764, as obtained from the DP estimate |^5] v\\/u± = 1.764(7). 
The DP estimate is derived from the scaling relation f\\/u± = 2/z using the earlier 
cited estimate for z. As can be seen the data for the dimer-trimer model is again fully 
compatible with DP critical behavior. 



3.3 Steady state behavior. 

Finally I have studied the steady state behavior of the density of sites. As mentioned 
earlier \p s ^ t — px\ should go to zero like a power-law at p\. In Fig. 5 I have plotted 
these quantities as a function of the distance from the critical point on a log-log scale, 
using p sat = 0.8935 for B particles and 0.1065 for isolated empty sites. The results 
were obtained by averaging over typically 100 independent samples. The number of 
time steps and system sizes varied from t = 1000, L = 64 far from pi to t — 500, 000, 
L = 512 closest to p\. The slope of the straight line is 0.59 as obtained from the DP 
estimate. From this figure it is clear that all densities have the same asymptotic power- 
law behavior and once again confirm that the model belongs to the DP universality 
class. 
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4 Summary and discussion 



Previous studies of various two-component models in two dimensions exhibiting a 

yielded 



continuous phase transition into a non-unique absorbing state \\L7\ , 
non-DP behavior. However, with the results reported in the previous section and in 
an earlier article | 23fl , it is now clear that of these both the dimer-trimer model [|17| 
and the CO-NO model [|IIJ belong to the DP universality class. Likewise, results 
for the one-dimensional versions of the pair contact process |2(| |2l| , dimer reaction 
model [[Hp and the threshold transfer process |22| clearly placed these models in the 
DP universality class. In all of these models the absorbing states can be uniquely 
characterized by the vanishing of a single quantity. More and more evidence thus 
confirm that the DP conjecture can be extended to such models, as first suggested in 
Ref. [21]. The exponent estimates for the various versions of the dimer-dimer model 
are not consistent with directed percolation. However, the typical value for /3 ~ 1/2 
(with an uncertainty of ~ 10%) is not very far from the DP value. In view of this, 
and the similarity of these models to the dimer-trimer and CO-NO models, it does not 
seem unlikely that the DD models also belong to the DP universality class. 
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Figure Captions 



Figure 1 Local slopes —S(t) (upper panel), rj(t) (middle panel), and z(t) (lower panel), 
as denned in Eq. |5| with m = 5. Each panel contains four curves with, from bottom to 
top, p = 0.3418, 0.3420, 0.3422 and 0.3424. 

Figure 2 The concentration of sites log 2 |/>x*(Pi> L)—px(pi, L)\ as a function of log 2 L. 
The slope of the straight line is P/u± = 0.81. Some of the densities have been scaled. 

Figure 3 The fluctuations in the concentration of active sites log 2 [x(pi> L)] vs log 2 L. 
The slope of the straight line is 7/z/_i_ = 0.39. 

Figure 4 The time before half the samples enter an absorbing state log 2 [r^(px, L)\ vs 
log 2 L. The slope of the straight line is v\\/v±_ = 1.764. 

Figure 5 Log-log plot of \px t (pi, L) — px\ as a function of the distance, \p — pi\, from 
the critical point pi = 0.3422. The slope of the straight line is /3 — 0.59. Some of the 
densities have been scaled. 
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